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Abstract 

> 
O 

CN New data imply that the average velocity of radio pulsars is large [1]. 

^sO Under the assumption that these data imply that a pulsar is born with an 

o 

"intrinsic" kick, we investigate whether such kicks can be a consequence of 
l^-) asymmetrical stellar collapse and explosion. We calculate the gravitational 

On 

wave (GW) signature of such asymmetries due to anisotropic neutrino radia- 
^ tion and mass motions. We predict that any recoils imparted to the neutron 

^ star at birth will result in a gravitational wave strain, hJ^ T , that does not go to 

a 

zero with time. Hence, there may be "memory" [2] in the gravitational wave- 
form from a protoneutron star that is correlated with its recoil and neutrino 
emissions. 

PACS numbers: 97.60.Bw, 97.60.Gb, 97.60.Jd, 95.30.Lz, 04.30.-w, 97.10.Wn 
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I. INTRODUCTION 



Recent data on pulsar proper motions [1], a recalibration of the pulsar distance scale [3], 
and a general recognition that previous pulsar surveys were biased towards low speeds [4] 
imply that many pulsars have high proper motions. Mean three-dimensional galactic speeds 
of 450±90 km s -1 have been estimated [4], with measured transverse speeds of individual 
pulsars ranging from zero to as much as ~1500 km s -1 . It has long been known that 
impulsive mass loss in a spherical supernova explosion that occurs in a binary can impart 
to the nascent neutron star a substantial proper motion that reflects its progenitor's orbital 
speed [5,6]. However, theoretical studies of binary evolution through to the supernova phase 
fail to reproduce velocity distributions with the requisite mean and dispersion [7,8]. 

Consistent evolutionary scenarios for the binary pulsars PSR BI534+I2 and PSR 
BI9I3+I6 [9], the observed misalignment of PSR BI9I3+I6's spin axis [10], a model for 
Her X-l [11], and the large eccentricities of Be/neutron star binaries [12] all buttress the 
conclusion that neutron stars are given an extra kick at birth. Recently, Caraveo [13] and 
Frail, Goss & Whiteoak [14] have identified pulsars with young supernova remnants (SNR's) 
in a majority of SNR's with putative ages less than 20,000 years, but only on the assumption 
that the pulsars have relatively high transverse velocities with an average value of ~500 km 
s" 1 . 

Related to the emerging SNR/pulsar associations are the new data on the "systemic" 
velocities of young supernova remnants [15,16]. It is observed that the "center-of-mass" 
velocity of oxygen clumps in four of these explosions is different from that of the local ISM 
by on average more than 500 km s -1 [16]. These data imply that the explosions themselves 
were asymmetric, but inhomogeneities in the ISM and a high progenitor speed can not yet be 
ruled out. It would be difficult to explain the jagged optical and IR line profiles of SN1987A, 
the intrinsic polarization of spectral features in SN1993J and SN1987A, and the oblateness 
of recent HST images of SN1987A without recourse to asymmetries during the explosion 
itself. In addition, Utrobin et al. [17] interpret the "Bochum" event in Ha in SN1987A with 
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a 10 3 M shard of 56 Ni moving at ~4700 km s 1 . Explosion asymmetries and the systemic 
SNR speeds may be (anti)correlated with the recoil of a nascent neutron star. 

II. RECOILS IN THE CONTEXT OF ASYMMETRIC COLLAPSE AND 

EXPLOSION 

In the past, an off-center (and rotating) magnetic dipole [18] and anisotropic neutrino 
radiation [19,21,22] have been invoked to accelerate neutron stars. It is the contention of this 
paper that asymmetries in the collapse and before and after the reignition of the supernova 
have the potential to impart to the core large recoils. That core-collapse supernovae are 
subject to Rayleigh- Taylor instabilities has been independently and convincingly demon- 
strated by at least three theoretical groups [23-25]. The pre-explosion core experiences a 
convective "boiling" phase behind the temporarily stalled shock, and the explosion, when 
it occurs, erupts in bubbles, fingers, and plumes. To date, these radiation/hydrodynamic 
calculations have been performed in only two dimensions (with axial symmetry) and with a 
variety of simplifying approximations. Nevertheless, this class of 2-D simulations can help 
theorists explore the potential role of hydrodynamics and neutrinos in imparting proper 
motions. The 2-D calculation of Burrows, Hayes, and Fryxell [23] was done assuming that 
the collapsing Chandrasekhar core and implosion were spherical. Even in this case, the core 
received stochastic impulses sufficient to shake and rotate the residue. 

Interestingly, recent hydrodynamic calculations [26] of convection during shell oxygen 
and silicon burning in massive stars and recent theoretical arguments [27] suggest that just 
before collapse asymmetries in density, velocity, and composition can be large. Furthermore, 
rotation might interact with convection to further distort the core [28]. The upshot might 
be an asymmetrical, aspherical collapse. If the amplitudes of the asymmetries in density or 
velocity are not negligible and if a significant low-order mode (l = 1) exists at the onset of 
collapse, the young neutron star can receive a large impulse during the explosion in which 
it is born. 
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To investigate this hypothesis, we conducted an exploratory and experimental calculation 
of aspherical collapse, using the code described in Burrows, Hayes, & Fryxell [23]. In this 
simulation, we artificially decreased by 15% the density of the Chandrasekhar core exterior 
to 0.9 M and within 20° of the pole. Such a perturbation in the 15 M progenitor that we 
used [29] amounts to a core- wide mass dipole anisotropy of less than 0.1% (and on any given 
mass shell of less than 2%). The initial total quadrupole anisotropy is 0.008 (normalized to 
2/3). This calculation was done in 2-D with azimuthal symmetry, but 6 ranged between 0° 
and 180°. Hence, the entire core, not just a wedge, was followed. To avoid severe Courant 
limitations, matter interior to 15 kilometers was followed in 1-D (radial). A 15% decrease in 
density is a bit larger than yet seen in the calculations of Bazan & Arnett [26], but we imposed 
no initial aspherical perturbation in velocity, despite the up to Mach 0.25 asymmetries they 
have seen. (Other calculations we have performed in which the initial asymmetry is solely 
in velocity, but of a similar magnitude, yielded quantitatively similar effects.) An essential 
point is that initial asphericities in this simulation grew during collapse, so that the mass 
column depths in various angular directions diverged. The matter collapsed at different rates 
in different directions, though pressure forces were transmitted in the angular directions as 
well that partially smoothed the deviations. The bounce was delayed on the side of the 
perturbation wedge and the resulting shock bowed out in the wedge direction. The accretion 
rates through this shock were highly aspherical. To avoid squandering computer cycles in 
what was merely a "proof-of-principle" calculation, we artifically hardened the emergent 
neutrino spectrum to facilitate an early explosion. The electron-type neutrino spectra were 
given an effective degeneracy factor,"??", of 3, in the upper range of the 1.5-3.0 normally 
encountered in fits to more realistic spectra [30]. Since neutrino heating drives supernovae, 
this ignited the explosion within 10 milliseconds of bounce. The subsequent explosion was 
aspherical not only due to the normal instabilities, but also due to the asphericity of the 
matter into which the explosion emerged and/or was driven. 

Figure 1 depicts the flow late in the explosion. The explosion erupted preferentially 
through the path of least resistance, i.e. in the direction of the wedge that we had imposed. 
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The wedge collapsed more slowly than the rest of the core. Since neutrino heating drives 
the explosion, matter heated near the neutrinosphere expands out as if from a reaction 
chamber. The protoneutron star residue receives a significant impulse a la the rocket effect. 
In this calculation, the center of the core is held fixed and we infer its recoil speed from the 
momentum it absorbs. The final core mass is ~1.2 M and its recoil speed is depicted in 
Figure 2. (The quasi-sinusoidal oscillation before bounce at 0.215 seconds seen in Figure 2 is 
the initial sound wave generated by the perturbation. Had the core not been anchored at the 
center, the anisotropy at bounce might have been even larger.) Figure 2 shows that the recoil 
speed grew almost monotonically and reached ~530 km s _1 . This is large, but only ~2% of 
the speed of the supernova ejecta. A small asymmetry in collapse clearly translated into an 
intrinsic kick, though the initial total momentum was almost zero. Just after bounce, the 
recoil is in the direction of the wedge, since the rest of the matter bounced first. Afterwards, 
as the shocked matter starts to squirt through the region of least resistance, the recoil 
changes sign and grows inexorably to its asymptotic value. The mass motions dominate 
the recoil. The contribution of neutrino radiation asymmetry to this kick is also depicted 
in Figure 2 and amounts to a total of ~I6%. The 77^, z/ T , and 77 T 's contribute 50% 
and the z/ e 's contribute 30% to the neutrino recoil. The magnitude of the dipole anisotropy 
of the radiation field varies from 1% to 10% (with an average of ~3% during the first 50 
milliseconds) and can change sign. Unfortunately, its behavior is not simple and at the end 
of the calculation, the dipole is +2%. On balance, it has the same sign as the term due to 
mass motions, since the neutrino fluxes are larger on the thin, or low-column, side. However, 
we have yet to demonstrate that this will be true universally. 



III. GRAVITATIONAL RADIATION SIGNATURE 

The impulse delivered to the core depends upon the dipole moments of the angular 
distribution of both the envelope momentum and the neutrino luminosity. The gravita- 
tional waveform depends upon the corresponding quadrupole moments. Using the standard 
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quadrupole formula, we derive for the neutrinos the expression, 

AG /■* 

Dh T J(t) = — / a(t')L v (t')dt', 

C J — oo 

where D is the distance to the supernova, h T z J is the transverse-traceless and dimensionless 
metric strain (z is the symmetry axis and the signal's directional dependence has been 
dropped), ait) is the instantaneous quadrupole anisotropy, and L v it) is the total neutrino 
luminosity [31,32]. Curiously, if a(t) is a constant, h T z J due to the neutrinos is proportional 
to the integrated neutrino energy loss! A burst of neutrinos translates into a rapid rise in 
h^J . Unfortunately, a(t) in the simulation was not constant, changed sign more than once, 
and on average was in the opposite direction to the neutrino's dipole term. Furthermore, 
the net contribution to h T 7 J of the neutrinos and the matter were of opposite signs. The 
magnitude of ait) averaged —0.02 during the first 50 milliseconds after bounce, but achieved 
—0.08 during the early shock break-out burst. The matter's contribution to both the dipole 
(recoil) and the quadrupole (GW) are of the same sign (as one might expect of an emerging 
"blob"). Figure 3 depicts the evolution of the strain(s) versus time. Due to the intense 
and anisotropic early neutrino burst, the neutrino contribution to changes in h T 7 J dominates 
during the first 20 milliseconds. Thereafter, that due to the expanding ejecta asserts itself 
(though at the end of this calculation, the neutrinos are continuing to contribute). This is 
true despite the fact that the neutrinos contribute only 16% to the recoil and is a consequence 
of the relativistic nature of the neutrinos. However, due to the more rapidly changing mass 
quadrupole around the time of anisotropic bounce, the total energy radiated (~ 1.1 X 10~ 9 
M c 2 ) is predominantly due to the mass motions (~30% is due to the neutrinos). 

The gravitational waves are radiated between 10 and 500 Hz. The "memory" in h T 7 J 
seen in Figure 3 is radiated at 10-100 Hz, just at the edge of the LIGO sensitivity range 
[33]. This memory is a distinctive characteristic of asymmetric collapse and explosion. Using 
reference [33] and Figure 3, we conclude that the 2'nd-generation LIGO will be able to detect 
the signal in Figure 3 at a signal-to-noise ratio of 10 from a core collapse anywhere in our 
galaxy. 
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A comparison with the results of Monchmeyer et al. [34] for rotating collapse is instruc- 
tive. In all their models, the asymptotic h T 7 J was zero and there was no memory. The 
maximum energy radiated in gravitational waves was for their Model A and was ~70 times 
that for our simulation. However, the maximum h T z J in their Model A was only ~3 times 
that depicted in Figure 3. Interestingly, their Model B experienced a centrifugal bounce at 
sub-nuclear density, radiated the same total wave energy as our model, but achieved only 
70% of its peak h^J . We see that even without rotation, an asymmetric collapse can result 
in appreciable gravitational wave emission. 

IV. CONCLUSIONS 

The major conclusion of this exploratory, but artificial, simulation is that a small initial 
core asymmetry can translate after collapse into an appreciable neutron star recoil. Other 
calculations that we have performed imply that this asymmetry /recoil correlation is generic. 
However, whether such asymmetries are themselves generic has yet to be demonstrated. 
This recoil is correlated with an appreciable gravitational wave signal with "memory," and 
both are correlated in a complicated fashion with the neutrino emissions. The simultaneous 
detection of these signatures would provide direct and unique information concerning the dy- 
namics of the supernova mechanism. Since accretion-induced collapse (AIC) is not preceded 
by the convective burning stages characteristic of the final hours of the core of a massive 
star, the initial asymmetries in the two contexts may be quite different. Consequently, the 
proper motions of AIC neutron stars and those of neutron stars from massive stars may be 
systematically different, with those of the latter being on average higher. This may be of 
relevance in the study of globular cluster neutron star binaries. 

The simulation highlighted in this paper is but an idealized study of the character of 
anisotropic collapse and explosion. It is in no way definitive. To explore these phenomena 
more credibly will require a realistically aspherical initial core, 3-D simulations, and much 
better neutrino transfer. Nevertheless, we have reproduced the correct magnitude of the 
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observed pulsar proper motions and have discovered potentially rich correlations between 
recoils, gravitational radiation, and neutrino emissions. Given the stochastic nature of the 
processes we highlight here, we expect that Nature provides not a single high kick speed, 
but a broad distribution of speeds. These will depend upon the degree and character of the 
initial asymmetries, the initial rotation structure, the duration of the delay to explosion, the 
progenitor density profiles, and chance. 
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FIGURES 




FIG. 1. A grey-scale rendering of the entropy distribution at the end of the simulation, about 
50 milliseconds into the explosion. Note the pronounced pole-to-pole asymmetry in the ejecta 
and the velocity field (as depicted with the velocity vectors). The physical scale is 2000 km from 
the center to the edge. Darker color indicates lower entropy and 9 = on the bulge side of the 
symmetry axis. 
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FIG. 2. The inferred recoil speed (in km s _1 ) imparted to the core versus time (in seconds) for 
the simulation highlighted in this paper. The initial momentum is approximately zero, but grows 
systematically after bounce in the direction opposite to the artificial wedge, cut into the core to 
mimic an asymmetry just before collapse. Shown are the total recoil (solid) and the contributions 
due to the neutrino emission anisotropy (dashed) and the ejecta motions (dotted). 
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FIG. 3. The gravitational wave strain, hJJ , times the distance to the supernova, D, versus 
time (in seconds). Core bounce is at 0.215 seconds. The total, matter, and neutrino waveforms 
are rendered with the solid, dotted, and dashed lines, respectively. 
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